arXiv:1502.07985vl [physics.optics] 27 Feb 2015 


The Effect of Experimental Parameters on Optimal Transmission of Light 

Through Opaque Media 

Benjamin Anderson 1 , Ray Gunawidjaja 1 , and Hergen Eilers* 1 

1 Applied Sciences Laboratory, Institute of Shock Physics, 

Washington State University, Spokane, WA 99210-169^\ 

(Dated: March 2, 2015) 

Spatial light modulator (SLM) controlled transmission of light through opaque media is a relatively 
new experimental method with wide applications in various fields. While there has been a surge 
in research into the technique there has been little work reported considering the effects of various 
experimental parameters on the efficiency of optimization. In this study we explore the effects of 
various experimental conditions on optimization and find that the intensity enhancement depends 
on the number of modulated channels, number of phase steps, feedback integration radius, beam 
spot size, and active SLM area. We also develop a model, based on the propagation of a Gaussian 
beam with a random phase front, to account for most of the measured effects. 
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I. INTRODUCTION 

Scattering materials (e.g. paper, paint, clouds, 
etc.) are generally viewed as hindrances to optical 
propagation. Light propagating through such a ma¬ 
terials behaves in a diffusive manner with scattering 
resulting in the amplitude and phase patterns be¬ 
ing destroyed [l|. However, in 1990 Freund theoret¬ 
ically showed that if you correctly shape the wave- 
front of the beam incident on the scattering system, 
the opaque system can be used as a lens or other 
high precision optical device |2]. 

The first experimental realization of an opaque 
lens was demonstrated by Vellekoop and Mosk [3j. 
Their technique involved using a liquid crystal on 
silicon (LCOS) spatial light modulator (SLM) to 
change the phase of the incident wavefront such 
that the modulated beam matches the transmission 
eigenmodes of the scattering sample. This leads to 
the light being focused through the sample onto a 
target area [3. 4] . This method of wavefront con¬ 
trol - using an LCOS SLM - has promising applica¬ 
tions for astronomical and biological imaging [5, 6], 
flourescene microscopy [3, 0, sub-diffraction limit 
focusing 0, fiol focusing and compression of ultra- 
short pulses [ill, m > spectral filtering HM3 , and 
light polarization control Ham. 

Our interest in the technique of SLM controlled 
optimal transmission is as a mechanism for verifying 
physically unclonable functions (PUFs) 0- PUFs 
are materials with a large number of random degrees 
of freedom that are practically impossible to recreate 
due to their inherent randomness pj)]. An example 
of a PUF is a scattering system, such as a nanopar- 
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tide (NP) doped polymer. In such a PUF the posi¬ 
tions of all the NPs are the degrees of freedom and 
since there are trillions of randomly distributed par¬ 
ticles in the scattering volume it is impossible to 
recreate the PUF; such that the scattering signature 
is identical. While the scattering signature (speckle 
pattern) is one way of authenticating/characterizing 
a scattering PUF, another possibility is to use SLM 
controlled transmission. For this application the 
SLM can be used to create a phase profile (challenge) 
that produces a specific transmission profile (key). 
If someone tampers with the PUF the coupling be¬ 
tween the phase profile and transmission profile will 
be broken, giving evidence of tampering. 

In order to use a SLM controlled transmission sys¬ 
tem for secure authentication of scattering PUFs we 
must first understand how different experimental pa¬ 
rameters affect the system. We therefore build a 
SLM controlled transmission setup and characterize 
- both theoretically and experimentally - the sys¬ 
tem’s optimization dependence on five different sys¬ 
tem parameters: SLM bin size, 6, number of SLM 
phase steps, M, active SLM area, L 2 , detector inte¬ 
gration radius, r, and the on-sample beam spot size, 
w. 

For comparing the influence of the different vari¬ 
ables on optimization we calculate/measure the in¬ 
tensity enhancement, 77, which is defined as: 


where / is the average intensity in the target spot 
after optimization and (Jo) is the ensemble aver¬ 
aged intensity in the target before optimization [ 3 ]. 
Figure |T] shows an example of the intensity pro¬ 
file both before-optimization and after-optimization. 
The pattern before optimization is a random speckle 
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pattern, while afterwards the pattern is a focused 
spot. By calculating the average intensity in the 
spot both before and after optimization we calcu¬ 
late the enhancement using Equation [lj 

The intensity enhancement due to SLM phase 
modulation has previously been modeled using an 
analog to electron conduction mil, which pre¬ 
dicts that the enhancement depends on the number 
of modulated SLM channels, TV, as 13]: 


V = ^(N- 1 ) + 1 . ( 2 ) 


From Equation [2] we see that the conduction model 
predicts that the enhancement should depnd only 
on the number of modulated SLM channels, with 
other system/sample parameters having no influence 
nil- In practice, however, experiments on a vari¬ 
ety of systems are found to give drastically different 
enhancements for similar N values [3|, [l7|, HH, l25l 1 2 81 ] . 
To account for these variations Yilmaz and cowork¬ 
ers developed a model to include detector noise into 
the optimization scheme [29]. Taking noise into ac¬ 
count, Equation [2] transforms into 
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the sample, with the beam waist being located at 
the incident surface. The incident electric field is 
therefore 


Ei(x,y) = E 0 e~( x2+y2 V° 2 ° (4) 

where Eo is the incident field strength and <jq is the 
beam’s Gaussian width. To model scattering, we let 
the beam width increase, a o cr, and introduce a 
random phase profile, With these transfor¬ 

mations the field exiting the sample is given by: 


E(x, y) = (5) 

We note here that this model of scattering is a sim¬ 
plistic approximation and does not reflect light prop¬ 
agation in real scattering media. However, due to 
it’s simplistic nature it allows us to perform a wide 
range of calculations, which are otherwise unfeasible. 

Assuming that the distance from sample to de¬ 
tector, Z, is much greater than the beam width, 
c j « Z , we can use Fraunhoffer diffraction the¬ 
ory to determine the beamprofile at the detector. In 
Fraunhoffer diffraction theory the diffracted electric 
field, Ed(x,y) [30j, is the Fourier Transform of the 
initial electric field: 


where R is the signal-to-noise ratio of the system. 

From Yilmaz and coworkers model, the enhance¬ 
ment should only depend on the number of SLM 
channels and the signal-to-noise ratio of the sys¬ 
tem. However, in our current work we find that the 
enhancement depends on the five system variables 
mentioned earlier, as well as the samples themselves. 
Given that our experimental results deviate from 
Equations [2] and [3] we propose a beam propagation 
model to account for the measured enhancement’s 
dependence on system parameters. The model- 
which we label the random phase Gaussian beam 
model (RPGBM)- is based on the propagation of a 
beam with a spatially random phase distribution and 
a Gaussian amplitude distribution. From both the 
model and experiment we determine phenomenolog¬ 
ical equations to describe the enhancement’s depen¬ 
dence on different system variables. 


II. MODEL 
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where k is the wavenumber given by k = 27r/A, with 
A being the wavelength of light. From the diffracted 
electric field we calculate the beam profile given by 


Id(x',y') = \E d (x',y')\ 2 . (7) 

We model the optimization process by introduc¬ 
ing a phase shift, ip(x,y), to the beam in the sam¬ 
ple plane, such that !p(x,y) represents the influence 
of SLM phase modulation. By systematically vary¬ 
ing !p(x,y) according to an optimization algorithm 
we can reproduce the same process used experimen¬ 
tally. With the addition of the SLM phase shift, the 
diffracted field becomes: 


A. Theory 

The random phase Gaussian beam model 
(RPGBM) treats the scattering sample as a “black 
box”, with the effect of scattering to introduce a ran¬ 
dom phase pattern to a Gaussian beam. We begin 
by assuming a TEMqo Gaussian beam incident on 


E, 


/ oo poo r 

/ dxdyE 0 exp l - 

-oo J — oo l 


(x 2 + y 2 ) 


, k 


-i®(x, y) + i-^[xx' + yy') + iip(x, y) 


( 8 ) 




3 



(a) Before (b) After 

FIG. 1: (Color Online) Camera image before (a) and after (b) optimization. The image before optimization 
is dim and random, while afterwards the beam is focused into a tight spot. The camera pixel size is 5.2 y m 
and the exposure time is 2 ms for the dim image and 0.127 ms for the optimized image. 


Using Equations [8] and [7] we can therefore model the 
detector feedback signal as we change i/;(x,y). 

While Equation [8] uses continuous Fourier Trans¬ 
forms, we use discrete Fourier transforms when 
performing computations. Discretizing Equation [8] 
gives 

N—l N — l 

Ed\n'm' = ^ ^ ^ ^ Eq GXp 
n=0 m =0 

2?r 

-i^n,m + ijjr(nn'+ mm , ) + itp ntm '>. (9) 

where we have substituted the x, y coordinates with 
integer values n, m such that: 


| ( n 2 + m 2 ) Ax‘ 

l 


Ax << cr, where a is the Gaussian width of the 
electric field. Using a random number generator 
each grid point is assigned a phase value between 
0 and 27r with the generated numbers having a uni¬ 
form probability distribution. The combination of 
the phase value, <F n?m , and Gaussian width, a de¬ 
fines the electric at the sample plane as 
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where ij; n,m comes from the SLM modulation and 
the peak field is defined as: 


x = nAx x' = rt Ax' 

y = rriAx y' = m Ax' 

with Ax being the grid spacing in the sample plane 
and Ax' is the grid spacing in the target plane given 
by 


Ax' = 


2ttZ 
NkAx’ 


( 10 ) 


Note that Equation [TO] implies an inverse relation¬ 
ship between distances in the sample and detector 
planes. 


B. Computational Details 

For our calculations we define a 1000 x 1000 grid 
with an isotropic grid spacing of Ax, such that 


E 0 = 


2Ax 2 


( 12 ) 


such that the total integrated intensity is unity. 

Using the sample-plane electric field from Equa¬ 
tion EH we model optimization as follows: 


1- V’ n,m is varied according to the chosen opti¬ 
mization algorithm and parameters. In this 
study we use a sequential bin-by-bin optimiza¬ 
tion method in which one bin is modulated at 
a time to find the optimal phase value, after 
which that phase is fixed for that bin. 


2. Given Ed ]n ',m' is calculated by taking 

the DFT of E n:Tn 


3. The intensity is found using Ed ]n ',m' and Equa- 
tion[7]and the average intensity, (/), in a target 
area of radius r is calculated. 
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4. A random number, 07 , which represents detec¬ 
tor noise, is added to the calculated average 
intensity with the random numbers having an 
average magnitude of (\cri\) = %J (/). 

5. The intensity with noise term, (I) + 07 , is then 
used as the feedback signal for the optimiza¬ 
tion algorithm. 

While we use a sequential bin-by-bin optimization 
algorithim in this study, the model’s optimization 
steps can easily be adapted for more complex algo¬ 
rithms, such as partitioning [4] and genetic optimiza¬ 
tion Hi Hi. 


C. Model Results 


1 . Bin Size 


We first model the effect of bin size, 6 , on op¬ 
timization by using an integration radius of 2 Ax' 
and three different numbers of phase steps, M = 
{2, 3, 20}. The calculations are performed both with 
and without noise with the result that the enhance¬ 
ment as a function of bin size depends on whether 
noise is present or not. Figure [2] shows a comparison 
of the enhancement as a function of inverse squared 
bin size for a calculation with and without noise. 
Without noise, the intensity enhancement is found 
to be proportional to a power function ( b~ 2 ) p where 
p < 1. However, when including noise in the calcula¬ 
tion we find that the intensity enhancement follows 
an exponential function, 


77 = 1 + 770 exp 



(13) 


where 1 + 770 is the asymptotic enhancement and 
abo is a shape factor, with a being a factor related 
to the active area of the SLM. Since the results of 
our calculations depend on whether or not noise is 
included, the remainder of our calculations we will 
include the effect of noise; as we can not completely 
eliminate noise experimentally. 

We next consider how changing the number of 
phase steps, M, influences the parameters in Equa¬ 
tion [13l Figure [3] compares the enhancement for dif¬ 
ferent number of phase steps, which we fit to Equa¬ 
tion m From our fits we find that the asymptotic 
enhancement increases as the number of phase steps 
increases, while the shape parameter remains con¬ 
stant. 



-2 -2 

Squared Inverse Bin Size b (units of Ax ) 


FIG. 2 : (Color Online) Modeled intensity 
enhancement as a function of inverse squared bin 
spacing. Without noise the enhancement follows a 
power function, while with noise the enhancement 
behaves as an exponential. 



FIG. 3: (Color Online) Modeled intensity 
enhancement as a function of inverse squared bin 
spacing for different M values. The enhancement 
follows an exponential function with the amplitude 
changing with M while the shape parameter 
remains constant. 





5 


2. Active SLM Area 

In the previous section we calculated the effect of 
changing bin size on optimization. This represents 
the first parameter which determines the total num¬ 
ber of controllable channels. The other parameter 
responsible for the total number of bins is the active 
SLM area, L 2 , with L being the active SLM side 
length. For modeling the effect of changing the ac¬ 
tive SLM area we use M — 10 phase steps, a bin 
size of b = 1 Ax, and three different radii: 1 Ax', 2 
Ax', and 5 Ax'. We calculate the enhancement as 
a function of active side length, shown in Figure 0J 
and find that it behaves as a Gaussian function: 


rj = 1 + rjo 


1 — exp 



(14) 


where 1 + 770 is the asymptotic enhancement and 
/3AL is the Gaussian width, with /3 being related 
to the bin size. Fitting the curves in Figure |4] we 
find that as the integration radius increases both the 
asymptotic enhancement and the Gaussian width 
decrease. This implies that to optimize a small ra¬ 
dius on the detector requires a much larger portion 
of the SLM to be active than in order to optimize 
a large target radius; which is expected given the 
inverse relationship between distances in the sample 
and detector planes. 


3. Phase Steps 

In addition to being able to change the number 
of controllable channels on the SLM, we also can 
vary the phase resolution of each channel. To model 
the effect of the number of phase steps on optimiza¬ 
tion we use an integration radius of r = 5 Ax', and 
four total bin numbers: N = {100,400,625,2500}. 
Figure [5] shows the enhancement as a function of 
phase steps, which is found to quickly saturate at 
around 10 phase steps for each N value. Qualita¬ 
tively this dependence can be understood because as 
the number of phase steps increases the phase reso¬ 
lution increases and at a certain point there will be 
diminishing returns in trying to attain higher phase 
resolutions. 

To understand this behavior we recall that the 
speckle pattern is an interference effect with the op¬ 
timization process attempting to match the phases 
of different beam portions to constructively interfere 
UM- This implies that the intensity in the target 
spot depends on a sum of interference terms of the 
form: 



FIG. 4: (Color Online) Modeled intensity 
enhancement as a function of active SLM side 
length. The enhancement is found to follow a 
Gaussian function. 


H n cos(<F n - V’n), (15) 


where A n is an amplitude factor, <F n is the phase 
before modulation, and i/j n is the contribution of the 
SLM given by 


Vvi 


2?r Q n 
M 5 


(16) 


with q n being an integer corresponding to the phase 
value giving the largest enhancement. While the ex¬ 
act functionality of the intensity is a complex sum 
over many such terms, we find that we can fit the 
enhancement as a function of phase steps using only 
one term giving a fit function of the form: 


r] = 1 + 770 cos 



(17) 


where 1 + 770 is the asymptotic enhancement and 0 0 , 
Ac/) are parameters which determine the shape of the 
function. 


4- Target Radius 

Thus far we have only considered how chang¬ 
ing the SLM properties affects optimization. How- 
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FIG. 5: (Color Online) Intensity enhancement as a 
function of the number of phase steps. 


ever, we also can control the detector’s parame¬ 
ters; most importantly, we can change the target 
integration radius. In order to model the effect of 
the target radius on enhancement we use M = 10 
phase steps and four different bin sizes such that 
N = {1000,2000,5000,10000}. Figure [6] shows the 
modeled enhancement which follows a double expo¬ 
nential as a function of the squared integration ra¬ 
dius: 


rj = 1 + Aie" r2/<T i + A 2 e~ r2/ ^ (18) 


where A\,A 2 are amplitude factors, and <7i, cr 2 are 
Gaussian widths. 

While the decrease in enhancement with increas¬ 
ing target area is expected, the functional form is 
surprising. To demonstrate this, we derive the ex¬ 
pected functional form by recalling that there is a 
finite amount of power, Po, that can be focused into 
the integration area. Assuming perfect enhancement 

- in which all the power is focused into the target 

- we would expect the maximum enhancement for a 
given radius to be 


Vmax ij ') 


(/o) 7 Tr 2 ’ 


(19) 


where (Jo) is the average intensity before enhance¬ 
ment. However, the RPGBM results are found to 
follow Equation [18] and not Equation [19] Currently 
the underlying physical principle determining this 
behavior is unknown and is an area of active re¬ 
search. 



2 2 

Integration Radius Squared r (units of Ax ) 


FIG. 6: (Color Online) Modeled intensity 
enhancement as a function of squared integration 
radius for four different bin numbers. The 
enhancement is found to follow a double 
exponential decay. 


5 . Beam Diameter 

The last system parameter we consider is the 
beam diameter. To model the effect of the beam 
diameter on the enhancement we use M — 10 phase 
steps, four different integration radii, and different 
bin sizes such that the beam diameter is always 10 
bins (i.e. b = 1 for a beam diameter of 10, b = 2 
for a diameter of 20, etc.). Figure [7] shows the en¬ 
hancement as a function of beam diameter for differ¬ 
ent integration radii. The enhancement is found to 
follow a peaked function where the peak location is 
dependent on the integration radius used. As the in¬ 
tegration radius decreases the beam diameter corre¬ 
sponding to peak enhancement is found to increase. 
This suggests an inverse relationship between the 
beam diameter and target spot size, which is consis¬ 
tent with the Fourier relationship between the two 
planes. 


D. Model Summary 

We model the process of SLM based transmis¬ 
sion optimization using a beam propagation model 
based on a Gaussian beam with a random phase pro¬ 
file. Using the model and a sequential bin-by-bin 
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FIG. 7: (Color Online) Modeled intensity 
enhancement as a function of beam diameter. The 
enhancement is found to be a peaked function, 
with the peak location being inversely related to 
the integration radius used. This result is 
consistent with the Fourier relationship between 
the sample and detector planes. 

optimization algorithim we optimize the diffracted 
pattern in a target area for varying systematic pa¬ 
rameters. We find that the optimization depends 
on all parameters tested which include: bin size, 
active SLM area, total number of phase steps, de¬ 
tector integration radius, and on-sample spot size. 
These results are different than those of previous 
models, which predict the enhancement to only de¬ 
pend on the number of modulated SLM channels 
(bins). @,[29]. 

In addition to our model predicting that the en¬ 
hancement depends on more parameters, it also pre¬ 
dicts a different dependence on the number bins 
used. To derive our systems dependnece on the 
number of bins we compare Equations [13] and [14] 
From these equations we find that the scale factors a 
and P relate the two equations with a = L/ A L and 
P = bo/b. Substututing the definitions of a and p 
into Equations [13] and [TT] we find that the model pre¬ 
dicts an intensity enhancement dependence on SLM 
bin size and active SLM area as: 

, = l+»exp|— (^l) }, (20) 

where bo is found to be independent of the number 
of phase steps and A L is found to decrease with in¬ 


creasing integration radius. Recalling that the num¬ 
ber of bins is given by N = (L/ 6 ) 2 , we can rewrite 
Equation [20] in terms of the number of bins, 

77 = 1 + 770 exp — T- 1 , ( 21 ) 

where Nq = (AL/bo) 2 - The bin number dependence 
in Equation [21] is drastically different than predicted 
by previous models 13. |29]. We hypothesize that 
this difference arises due to the inclusion of beam 
propagation effects and enhancement saturation due 
to detector noise. 

The other difference between the RPGBM and 
previous models, is that the RPGBM predicts that 
the number of SLM phase steps, detector integration 
radius, and the on-sample beam spot size also affect 
the intensity enhancement. The dependence on the 
number of phase steps arises due to optimization be¬ 
ing related to controlled interference, while the influ¬ 
ence of the integration radius and on-sample beam 
spot size occurs due to the diffractive nature of the 
enhancement phenomenon. Since the sample and 
detector planes are related via a Fourier Transform 
the effect of changing distances in one plane directly 
affects distances in the other plane. 

Finally, one of the major simplifications of 
RPGBM is using a smooth Gaussian amplitude for 
the scattered beam, whereas real scattering results 
in both the amplitude and phase of the beam being 
modulated. Experimentally we find that the ampli¬ 
tude of the beam leaving the scattering sample is 
approximated by a Gaussian with an additive noise 
term. When performing calculations we find that 
adding a noise term has negligible effect on the func¬ 
tional form of the optimization. This is due to the 
model having the SLM only affecting the phase of 
the scattered beam and not its amplitude. In reality, 
however, the scattered beam’s phase and amplitude 
are coupled due to the material and we expect that 
using an SLM to change the incident beam’s phase 
will result in a small change in the scattered am¬ 
plitude. To account for this effect we are currently 
extending the RPGBM with a model of scattering 
based on transmission eigenchannels and random- 
matrix theory [4, [ 2 (J [ 2 l|, 24]. 

III. EXPERIMENTAL METHOD 

We experimentally measure the effect of various 
system parameters on optimization using a con¬ 
trolled transmission optical setup, which consists of 
a high-speed LCOS SLM from Boulder Nonlinear 
Scientific, a Coherent Verdi V10 Nd:YV 04 laser, a 
high speed Thorlabs CMOS camera ( 8 -Bit,pixel size 
of 5.2 /im, variable exposure time from 37 fi s to 2 


















IV. RESULTS AND DISCUSSION 


90:10 BS 



FIG. 8 : (Color Online) Optical setup schematic. 


ms), and various focusing and polarization optics. 
Figure [ 8 ] shows a schematic of the system. 

The beam from the laser is sampled by a 90:10 
beamsplitter (BS) and then expanded by a factor of 
3.75 x and passed through a half-waveplate (HWP) 
polarizer pair to control beam intensity. The ex¬ 
panded beam is then reflected by a beamsplitter 
onto the SLM which modulates the phase of the 
laser beam. After modulation the beam is focused 
onto the sample using a 20 x high working distance 
(HWD) objective (WD = 20 mm), with the scat¬ 
tered light being collected onto a CMOS camera by 
a 5x HWD objective (WD = 37.5 mm). 

To keep optimization times managable we bin the 
SLM pixels into N bins with each bin having an 
edge size of b = L/7V, where L is the total number 
of active pixels on a side. The bins are optimized 
using a sequential optimization algorithm II in 
which each bin’s phase value is updated through M 
phase steps of size A <fi = 2tt/M. After each update 
the camera takes an image which is used to calcu¬ 
late the intensity within the target area. After all 
M steps are completed, the bin’s phase is fixed to 
the phase value corresponding to the largest inten¬ 
sity measured. This procedure repeats for all bins 
until an optimized phase pattern is displayed on the 
SLM. Using this algorithm with our setup we achieve 
iterative rates of 160-180 Hz. 


To characterize the controlled transmission setup 
we systematically vary the five system parameters: 
bin size, active SLM area, number of phase steps, 
target area, and beam spot size. We measure each 
dependence by varying one parameter, while holding 
all other parameters fixed, and measure the intensity 
enhancement. To obtain better statistics we perform 
ten optimization runs for each parameter set and 
find the average enhancement and error from the 
ten runs. 

In order to separate which effects are due to 
the systematic parameters and which are due 
to the opaque sample, we perform measurements 
on five different sample types: Zr 02 NP-doped 
polyurethane (PU), Zr 02 NP-doped polyepoxy (PE), 
Y 2 O 3 pressed ceramic, ground glass, and printer pa¬ 
per. From these measurements we find that the func¬ 
tional form of the enhancement as a function of sys¬ 
tem parameters is independent of sample type, with 
the different samples only affecting fit parameters 
(e.g. peak enhancement, shape parameters). This 
suggests that the measured dependencies are a func¬ 
tion of the optical setup and not the samples. Since 
the functional forms are consistent across samples; 
we present in the following sections the enhancement 
measured using ground glass, as it’s speckle pattern 
is found to be the most stable over time and it pro¬ 
duces the largest enhancements. 


A. Bin Size 

The first parameter we investigate is the bin size. 
We measure the intensity enhancement at 12 differ¬ 
ent bin sizes using a beam spot of diameter « 350 
/im, an integration radius of 2 px, and three different 
total number of phase steps, M = {8,10,16}. Figure 
[9] shows the enhancement as a function of squared 
inverse bin size, which is found to follow a function 
of the form 


7? = l+ % (l-e- b > 2 ), (22) 

where 1 + 770 is the asymptotic enhancement and bo 
determines the enhancement’s shape. We find the 
parameters 770 and bo for each M value by fitting 
the curves in Figure [9j which are tabulated in Table 
HI The asymptotic enhancement is found to increase 
with the number of phase steps, while the shape fac¬ 
tor, bo , is found to be constant within uncertainty. 
These results are functionally consistent with the 
RPGBM results. The difference in the magnitude of 
enhancement between model and experiment is due 
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FIG. 9: (Color Online) Intensity enhancement as a 
function of squared inverse bin size. The 
enhancement is found to follow an exponential 
function, which is consistent with the RPGBM 
results. 


M 

Vo 

bo 

8 

10 

16 

23.4 ±4.9 
50.1 ±7.2 
146 ± 14 

21.4 ±3.5 
19.9 ± 1.9 

19.5 ±1.3 


TABLE I: Fit parameters from Equation [22] for the 
intensity enhancement as a function of bin size. 


to three factors: 1) imperfect matching of parame¬ 
ters between experiment and modeling, 2) different 
noise levels, and 3) a divergence between the model 
and experiment related to the number of phase steps 
(which is discussed later in the paper). 


B. SLM Cropping 

Next we test the effect of changing the active SLM 
area using a spot diameter of ~ 200 /im, a bin size 
of 6 = 8 px, M = 8 phase steps, and three different 
integration radii, r = {1,2,5}. We first measure the 
enhancement with the full SLM active, after which 
we “shut off” the outer rows/columns - such that 
the active area is always a centered square - and 
perform optimization again. We continue to do this 
until only a 4 bin x 4 bin (32 px x 32 px) area re¬ 
mains active. Figure [TO] shows the intensity enhance¬ 
ment as a function of quartic active side length, L 4 . 
From the figure we find that the enhancement, as 


r 

Vo 

A L 

1 

2 

5 

32.10 ±0.80 
24.36 ± 0.25 
5.94 ±0.11 

232.6 ±3.6 

168.3 ± 1.9 

156.4 ±3.6 


TABLE II: Fit parameters from Equation [23] for 
the intensity enhancement as a function of active 
area. Both 770 and A L decrease with increasing 
integration radius, which is consistent with the 
RPGBM. 


a function of active length, L, is found to behave 
functionally as 


A L ) 

where A L is a width parameter and 1 ± 770 is the 
asymptotic enhancement. From Equation [23] we find 
that the explicit dependence of the enhancement 
on active side length is different than predicted by 
RPGBM (i.e. Gaussian in L for the RPGBM and 
Gaussian in L 2 for experiment). 

A possible explanation for this discrepancy is re¬ 
lated to how the RPGBM treats SLM cropping ver¬ 
sus the real world implementation. In the RPGBM, 
cropping of the active SLM area is implemented by 
shutting off modulation in grid points that repre¬ 
sent the sample’s exit plane. However, in reality 
SLM cropping shuts off modulation of portions of the 
light incident on the sample. The modulated light 
incident on the sample is then transmitted through 
the sample with various spatial components inter¬ 
fering. This interference is not accounted for in the 
RPGBM, which could lead to the divergence be¬ 
tween the RPGBM and experiment. 

Despite the explicit dependence on L being differ¬ 
ent, the general dependence - enhancement increas¬ 
ing to a constant value as L increases- is consistent. 
Additionally, we find that the behavior of the en¬ 
hancement’s dependence as the integration radius 
changes is also consistent. To demonstrate this con¬ 
sistency we fit Figure [lO] to Equation [23] and find 
770 and A L for each integration radius tested. Ta¬ 
ble [IT] compiles the fitting results. Both the ampli¬ 
tude, 770 , and width parameter, AT, are found to 
decrease with increasing integration radius, which is 
predicted by the RPGBM. 



77 = 1 ± 770 


1 — exp < 


C. Phase Steps 

The final SLM parameter we vary is the number 
of phase steps used during optimization. For these 
measurements we use a spot diameter of « 350 fi m, 
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FIG. 11: (Color Online) Intensity enhancement as 
a function of phase steps. 


FIG. 10: (Color Online) Measured intensity 
enhancement as a function of the quartic active 
side length, L 4 . The enhancement is found to 
follow an exponential function, which is different 
than predicted by the RPGBM. 

an integration radius of 2 px, and three bin sizes: 
b = 16 px, b = 32 px, and 6 = 64 px. We find that 
the enhancement depends on the number of phase 
steps, M, as: 

^l + ^co^yy, (24) 

where 1 ± r]o is the asymptotic enhancement and p 
is an exponent which controls the shape of the func¬ 
tion. While the RPGBM predicts that p = 1, from 
fitting we find that p > 1 and decreases as N in¬ 
creases, as shown in Table ini This result is unex¬ 
pected and the underlying mechanism is currently 
unknown. 

One possible explanation is that realistic samples 
may complicate the coupling between the modulated 
phase incident on the sample and the phase exiting 
the sample such that only a fraction of the light ex¬ 
iting the sample has an optimized phase. The result 
of having less control would be to decrease the en¬ 
hancement, which is consistent withp > 1. Also, we 
would expect effects due to the sample to decrease as 
N increases, since the size of the modulated area de¬ 
creases and gives a greater control over transmission 
through the sample. This increased control would 
cause p to decrease, which we see experimentally. 


N 

Vo 

P 

64 

256 

625 

1024 

65.5 ±6.4 
78.8 ± 6.4 
95.3 ±4.1 
116.0 ±8.2 

61.0 ±5.3 
59.2 ±6.5 
53.5 ±8.6 
29.0 ±3.5 


TABLE III: Fit parameters from Equation [24] for 
the intensity enhancement as a function of the 
number of phase steps. 


D. Target Radius/Area 

Thus far we have only considered the effects of 
SLM parameters. At this point we turn our at¬ 
tention to the enhancement’s dependence on tar¬ 
get radius. We measure the intensity enhance¬ 
ment using a beam spot size of 380 /rm, M = 
32 phase steps, four different numbers of bins, 
N = {256,625,1024,2025}, and nine integration 
radii/areas. The measured enhancement as a func¬ 
tion of integration radius, shown in Figure [12j is 
found to behave as the sum of two Gaussians. This 
behavior is identical to the RPGBM. 

We also determine amplitudes and Gaussian 
widths as a function of N by fitting the measured 
enhancement as a function of integration radius to 
Equation [T8l Table HVl compiles the fit results. The 
amplitudes are found to increase with bin number - 
consistent with our other results - and the widths 
are found to decrease as the bin number increases. 
Additionally, the Gaussian widths appear to reach a 
constant value as the number of bins increases, with 
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FIG. 12: (Color Online) Intensity enhancement as 
a function of integration radius measured using 
ground glass with 32 phase steps, a spot size of 380 
/i m and four different total number of bins. The 
enhancement is found to behave as a double 
exponential which is consistent with the RPGBM. 


N 

A! 

CFl 

a 2 

CT2 

256 

625 

1024 

2025 

22 ±13 
67.1 ±3.3 
129 ± 16 
190 ± 19 

27.1 ±5.5 

19.1 ±1.7 
15.85 ±1.66 
16.10 ±0.91 

56 ±12 
344 ± 17 
312 ± 33 
496 ± 71 

7.4 ± 1.5 
4.23 ±0.18 
3.95 ± 0.38 
3.77 ±0.35 


TABLE IV: Fit parameters from Equation [T8l for 
the intensity enhancement as a function of 
integration radius. 


the widths for N = 1024 being within uncertainty 
of those for N = 2025. 

As with the RPGBM model’s results, the under¬ 
lying mechanism behind the experimental enhance¬ 
ment’s target area dependence is currently unknown. 
Given that both the model and experiment have 
the same functional dependence, suggests that the 
physical phenomenon responsible is related to beam 
propagation effects. Currently we are performing 
further experiments and modeling in order to better 
understand how other parameters, besides the num¬ 
ber of bins, affects the enhancement’s dependence 
on integration radius. 



FIG. 13: (Color Online) Peak enhancement as a 
function of position along the optical axis, where 
z = 0 is the focal point of the focusing lens. The 
spot size at z = 0 is ~ 0.9 /rm. 


E. Spot Size 

The last system parameter we vary is the on- 
sample beam spot size. To measure the enhance¬ 
ment’s spot size dependence we use N = 1024 bins, 
M — 16 phase steps, and integration radii of 2 px, 20 

px, 30 px. We first measure the enhancement with 
the sample positioned within the focal length of the 
focusing objective such that the on-sample beam di¬ 
ameter is 600 fi m. After the initial measurement 
we systematically translate the sample and measure 
the enhancement at fixed z positions until the sam¬ 
ple translates through the focal point and reaches an 
on-sample beam diameter of 600 fi m again. Figure 
fT3l shows the peak enhancement as a function of po¬ 
sition along the optical axis, where z = 0 is the focal 
point of the focusing lens. From Figure [13] we find 
that the enhancement is symmetric about the focal 
point, with the peak value occurring at a nonzero 
distance from the focal point. 

We can convert the ^-position into the spot diam¬ 
eter, 2re, using ray matrix Gaussian beam propaga¬ 
tion. Assuming that the beam incident on the fo¬ 
cusing lens is at/near it’s waist, the Gaussian width, 

re, at position 2 is given by: 

W(z) = ^2—^4/2(/ + zf + (25) 
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FIG. 14: (Color Online) Measured intensity 
enhancement as a function of the on-sample spot 
size. The enhancement is found to follow a peaked 
function with the peak location being inversely 
dependent on the integration radius, consistent 
with the RPGBM result. 


where wo is the beam diameter at the focusing lens, 
/ is the focal length of the lens, and k = 2 tt/X where 
A is the wavelength of light. Note that Equation [25] 
is symmetric about the focusing lens’s focal point. 
Therefore we average the enhancement measured on 
both sides of the focal point to find the intensity 
enhancement as a function of spot diameter which 
is shown in Figure [lU 

The experimentally measured intensity enhance¬ 
ment is found to peak at a nonzero spot diameter, 
with the diameter corresponding to peak enhance¬ 
ment increasing as the integration radius decreases. 
Additionally the width of the peaked function is 
found to decrease as the integration radius decreases. 
These behaviors are consistent with the prediction 
of the RPGBM. 


F. Different Samples 

The experimental enhancement’s dependence on 
the different systematic parameters diverges from 
both previous models (Equation [2] and [3]) (3[l29j| and 
the RPGBM. While the RPGBM predicts similar 
behavior to experiment, some of the dependencies 
are functionally different. The most likely source 
of divergence between experiment and the RPGBM 


is the model’s treatment of the sample as a “black 
box”. This treatment predicts that the enhancement 
is independent of the sample properties (e.g. sample 
thickness and scattering length). 

While a precise characterization of the enhance¬ 
ment’s dependence on sample parameters is beyond 
the scope of this paper, we consider a simple di¬ 
rect comparison between six different samples: pa¬ 
per, ground glass, Y 2 O 3 ceramic, and three formu¬ 
lations of Zr 02 NP-doped polymers. We perform 
optimization using the same experimental parame¬ 
ters (b = 16 px, M = 32, r = 2 px, w = 250 fi m) 
at five different points on each sample to find the 
spatially averaged intensity enhancement, which is 
tabulated in Table [Vl From Table [V| we see a wide 
variation in enhancement when using different sam¬ 
ples, with the largest enhancement being 36 x larger 
than the smallest. 

We deduce several possible factors from these pre¬ 
liminary measurements which may affect the en¬ 
hancement: the sample persistence time and scatter¬ 
ing length. The first factor, the persistence time, is a 
measure of how long a sample will produce the same 
speckle pattern; which directly affects how well an 
SLM system can optimize transmission SI. From 
our measurements we find that the ground glass 
and NP samples have stable speckle patterns over 
a period of days, while the paper’s speckle pattern 
changes in tens of minutes. This results in paper 
having a relatively low enhancement despite being 
the thinest sample. While the persistence time is 
important to optimization, it does not explain the 
wide variation in the most stable samples. For these 
samples we see that the enhancement is largest for 
large scattering lengths (ground glass) and smaller 
for small scattering lengths (NP samples). 

The result linking the enhancement to the scat¬ 
tering length is important to the long-term goals 
of our study, which are to use optimal transmission 
as a method of authenticating NP-doped polymeric 
PUFs. With this in mind we note that for the two 
different 10 wt% Zr 02 NP-doped polymers the en¬ 
hancements are within uncertainty of each other, de¬ 
spite the samples having different host polymers and 
different thickness. Additionally, comparing 1 wt% 
and 10 wt% Zr 02 NP-doped PU, we find that the 1 
wt% enhancement is almost 10 x as large as for the 
10 wt%. These results suggest that the polymeric 
hosts have a negligible effect on the optimization 
process, while the particle concentration has a large 
inverse effect. This can be understood according to 
Mie’s scattering theory which predicts that the scat¬ 
tering length - for a system of scattering spheres - is 
inversely proportional to the scatterer density (3lj . 
This implies that as the concentration increases the 
number of scattering events increases. We hypothe- 
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Sample 

Thickness (/mi) 

Scattering Length (/im) 

Enhancement 

Ground Glass 

1564 ± 75 

970.7 ±2.1 

172 ± 12 

Paper 

85.1 ±6.4 

(2.654 ± 0.026) x 10“ 3 

9.4 ±1.0 

10 wt% Zr0 2 NP/PU 

867 ± 67 

4.11 ±0.28 

5.50 ±0.45 

10 wt% Zr0 2 NP/PE 

1036 ± 50 

3.9 ±1.2 

4.79 ± 0.49 

1 wt% Zr0 2 NP/PU 

959 ± 37 

50.8 ±3.1 

44.5 ± 1.2 

Y 2 0 3 Ceramic 

358 ± 44 

(2.944 ± 0.029) x 10“ 2 

12.8 ± 1.3 


TABLE V: Maximum enhancement obtained for different samples using system parameters of b = 16 px, 
M = 32, r = 2 px, w = 250 /im. There is a 36x difference between the smallest and largest enhancement. 


size that this increase in scattering events decreases 
our ability to control transmission of light through 
the system. While the exact mechanism for this 
effect is unknown, Mosk has proposed that by in¬ 
creasing the number of scattering events, either by 
increasing thickness or concentration, will result in 
noise playing a larger effect in optimization [32], re¬ 
sulting in the enhancement decreasing [29]. Further 
work is required to better understand this effect. 

V. CONCLUSION 

We systematically measure the dependence of op¬ 
timal transmission on five different system variables: 
SLM bin size, number of SLM phase steps, active 
SLM area, detector integration radius, and the on- 
sample beam spot size. From our measurements 
we find optimization to depend on all five system 
variables as well as the characteristics of the sample 
used for optimization. These results are contradic¬ 
tory to previous models of universal optimal trans¬ 
mission which proposed that optimization is only de¬ 
pendent on the number of SLM channels used |3|, |24| 
or the number of SLM channels and the signal-to- 
noise ratio [2Sj. To understand the nature of these 
contradictory results we develop a model based on 
the propagation of a Gaussian beam with a random 
phase profile. 

We find that the model is mostly consistent with 
experimental results, with the effects of beam prop¬ 
agation primarily arising due to the Fourier relation¬ 
ship between the sample and detector planes. This 
relationship leads to the beam diameter and active 
SLM area being inversely related to the target area 
to be optimized. This implies that to best optimize 
a small target radius a large beam spot size and 
large active SLM area are required, while to opti¬ 
mize a large target area requires less of the SLM to 
be active and a smaller spot size. We also find that 
including noise effects into the model leads to the 
correct enhancement dependence on the number of 
bins. 


While the RPGBM and experiment are mostly 
consistent, there are still some key variations: 
namely, the enhancement’s dependence on the num¬ 
ber of phase steps and active SLM area. These devi¬ 
ations most likely arise due to the RPGBM’s treat¬ 
ment of the sample and scattering, which assumes 
that the sample does not affect the optimization 
process and that the scattered phase and amplitude 
are independent. These assumptions, however, are 
incorrect as we see from experiment that the sam¬ 
ple does affect the optimization process (primarily 
due to the persistence time and the concentration 
of scatterers) and that the scattered amplitude and 
phase are correlated by scattering within the sam¬ 
ple. To address these shortcomings we currently are 
in the process of performing systematic studies of 
how different sample properties affect the enhance¬ 
ment as well as working on extending the RPGBM 
to more accuaretly represent scattering within the 
sample. To extend the RPGBM we are utilizing two 
different computational methods to describe scat¬ 
tering: random matrix transmission eigenchannels 
H, I 21 I, [24| and monte carlo scattering simulations 
33l-[36|. The scattering models will be implemented 
into the RPGBM as follows: for each phase mask 
tested the beam will be passed through the scatter¬ 
ing algorithim to determine the beam amplitude and 
phase exiting the sample. As both algorithims de¬ 
pend on many iterations using random numbers we 
will run the same phase mask multiple times to de¬ 
termine an ensemble average amplitude and phase. 
These ensemble averaged exit beams will then be 
transformed using Equation [8] to determine the in¬ 
tensity in the detector plane. This process will be 
repeated according to the optimization algorithim 
until the optimal phase mask is determined. 
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